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We report on the hadron mass spectrum obtained on a 16^ X 40 lattice in full QCD at /3 = 5.7 usins two flavors 
of staggered fermions with ma = 0.01. We study the effective mass plateaus for different sized sources. Our mass 
results are slightly lighter than our earlier 16^ x 32 calculation. The Landau gauge A is quite different from the 
Coulomb gauge A. 



1. INTRODUCTION 

The hadron mass spectrum is one of the fun- 
damental quantities to calculate in lattice QCD. 
The Columbia group has done hadron mass cal- 
culations with two flavors of staggered dynami- 
cal quarks for ma = 0.01,0.015,0.02,0.025 on a 
16^ X 32 volume and ma = 0.01,0.025 on 32^, 
both at P — 5.7 0-1). We also investigated 
ma = 0.004 on 16^ x 32 at ^ = 5.48 g. 

In this note, we report a new study on a IG'^ x 40 
lattice with ma = 0.01 at f3 — 5.7. The increased 
time dimension from Nt = 32 to Nt =40 gives 
longer plateaus in the effective mass plots so that 
hadron masses can be extracted more reliably. 
We systemically study the source size effects on 
the determination of the mass values. Two non- 
local A operators have also been added. They 
are measured in both Coulomb gauge and Lan- 
dau gauge. 

2. SIMULATION 

Table lists the parameters of our calcula- 
tion. We use the 'R'-algorithm of Gottlieb, et. al. 
j5| for our evolution along with their notation. 
We have collected 3100 microcanonical time units 
during a 5 month nm on the 256-node Columbia 
parallel computer. Hadron masses are measured 
every 6 time units. For a single gauge conflgura- 
tion, sources are set on 5 different time slices to 



perform an average over time slices (ACTS). 
Table 1 

Simulation parameters 



volume 

'^dynamical^ 

length 

trajectory length 
step size 

CG stopping condition 
minutes / trajectory 

hadron source types 

hadron source sizes 
valence quark masses 



16^ X 40 
5.7 
0.01 
250-3100r 

0.5r 
0.0078125T 
1.13 X 10-5 

12.5 

wall, 2Z 
wall, Z 
16^, 123,33,43, point 
0.004, 0.01, 0.015, 0.02 
0.025, 0.05, 0.07 



*This work was done in collaboration with Norman H. 
Christ, Robert D. Mawhinney, Igor Arsenin, Shailesh 
Chandrasekharan, Weonjong Lee and Decai Zhu and sup- 
ported in part by the Department of Energy. 



For hadron mass calculations, our 2Z source is 
set on a fixed t slice and is non-zero when only 
all {x,y,z) are even; our Z source is non-zero 
for all {x,y,z). Previously our 2Z and Z wall 
sources were non-zero across the entire spatial 
volume. In this run, however, we also set sources 
on smaller cubic spatial volumes. Besides the 
dynamical quark mass, we also calculate quark 
propagators using different valence quark masses 
so that we can make hadron propagators out of 
different sets of quark masses. Table ^ also lists 
all of our hadron source sizes and valence quark 
masses. Out of all the possible combinations of 
source types, sizes and valence quark masses, we 
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have a total of 19 combinations for mesons, 21 
combinations the nucleon and 10 combinations 
for the A. 

Figure |^ shows the evolution of (xx) during 
our run. Discarding the first 250 time units for 
thermalization gives (xx) — 0-0276(2). This is 
consistent with the value 0.0277(3) from our ear- 
lier 16^ X 32 run 



fact, in our earlier 32**, f3 — 5.7, ma = 0.025 calcu- 
lation with a 32'^ source, all the effective masses 
rose from below and the nucleon probably never 
reached a plateau. For our current data, all the 
effective masses from different sized sources con- 
verge to the same plateau. It is clear that 16^ 
is the optimal source size for this calculation be- 
cause it gives the longest plateau. 
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Figure 1. The evolution of the chiral Condensate. 



3. RESULTS 

We see quite significant source size effects on 
the effective mass plots for hadrons. Figure ^ 
and H show from top to bottom the effective 
mass plots for N, p, tt on 16^ and 8'^ sized sources. 
For the 16'^ source, the effective masses are quite 
flat as the distance varies. When the source 
size becomes smaller, the effective masses reach 
a plateau from above. The smaller the size of 
the source, the higher the value of the effective 
masses at short distances. This behavior is seen 
with all our different sources: 4^,8^,12^,16^. In 
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Figure 2. iV, p, tt effective mass 
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Tabic g is a list of hadron mass values for our 
current 16^^ x 40 calculation. Also listed for com- 
parison are our earlier 16^ x 32 results. All our 
fits are correlated fits with errors determined from 
the jackknife method. The pion mass is obtained 
from a single state fit with the 2Z source. All 
the other mesons and the nucleon come from two 
state fits with the 2Z source. Our A's come from 
two state fits with the Z source. Ai and Aq are 
two equivalent lattice operators corresponding to 
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Figure 3. N, p, n effective masS' 




5 .1 

Figure 4. Aq effective masses 
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the two A2 operators in ref ||], respectively. All 
our sinks are local except for the A's which are 
necessarily non-local in the staggered fermion for- 
malism. Figure ^ is an effective mass plot for the 
Coulomb gauge Aq with a 16'^ source. 

All our mass values are some what lighter than 
that found in the earlier 16^ x 32 run. The nu- 
cleon is especially light, about 3 standard devi- 
ations lower than the earlier result. This is not 
because of different fitting ranges; fitting our data 
with distances up to 16 instead of 20 gives essen- 
tially the same mass spectrum. The two runs 
have about the same statistics. The earlier one 
has more averaging on a single configuration while 
the current one has more independent configura- 
tions. Except for the possibility that we are see- 
ing a finite temperature effect (With antiperiodic 
boundary conditions in the time direction, our 
larger Nt lowers the temperature from 53 MeV 
to 44 MeV.), one may argue that in both of the 
runs, the errors are underestimated and that our 
simulation time is not sufficiently large to accu- 



rately compute the effects of autocorrelations. 

The A masses measured in Landau gauge are 
quite different from those measured with both the 
source and the sink time slices in Coulomb gauge. 
They are about 3 standard deviations lower. In 
fact, their masses are in closer agreement with 
the nucleon than the Coulomb gauge A's. This 
raises a question about the validity of using Lan- 
dau gauge, non-local operators in general, since 
the gauge-dependent correlation between differ- 
ent time slices in Landau gauge may alter the 
signal one wants to obtain. 

Setting the p mass to 770 McV to fix the scale, 
we find = 436(6) MeV, m^^ = 502(9) MeV, 
ruN = 1153(20) MeV and uia = 1290(40) MeV 
(averaging over our Coulomb gauge Aq and Ai). 
Our lattice spacing is = 1750(20) MeV and 
the lattice size is La = 1.80(2) fm. 
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Table 2 

Hadron mass results: (3 = 5.7, ma = 0.01. 



volume 


16^ X 32 


16'^ X 40 


length 


1050 2700T 


250 3100r 


measured every 


5t 


6t 


# of measurements 


330 


475 


AOTS 


32 


5 




0.0277(3) 


0.0276(2) 




0.252(3) 


0.249(2) 




0.295(3) 


0.288(4) 


rupa 


0.454(4) 


0.440(5) 


rup^a 


0.457(7) 


0.440(6) 


mNd 


0.692(6) 


0.659(9) 


rriN^a 


0.795(10) 


0.767(18) 


m^a 


0.441(14) 


0.435(10) 


mAiCi 


0.586(14) 


0.576(15) 


mBd 


0.610(9) 


0.609(17) 


rriN/mp 


1.527(11) 


1.499(13) 


mAoa(C) 




0.753(23) 


m-Aia(C) 




0.721(19) 


mAoa(L) 




0.688(15) 


mAia(L) 




0.659(22) 



The fits have tmin = 9 for the tt and tmin = 6 
for all other particles. The jdof for all fits are 
< 2. 
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